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Abstract 

We analyze here the issue of local versus the global visibility of a singularity that 
forms in gravitational collapse of a dust cloud, which has important implications for 
the weak and strong versions of the cosmic censorship hypothesis. We find conditions 
as to when a singularity will be only locally naked, rather than being globally visible, 
thus preseving the weak censorship hypothesis. The conditions for formation of a black 
hole or naked singularity in the Szekeres quasi-spherical collapse models are worked 
out. The causal behaviour of the singularity curve is studied by examining the outgoing 
radial null geodesies, and the final outcome of collapse is related to the nature of the 
regular initial data specified on an initial hypersurface from which the collapse evolves. 
An interesting feature that emerges is the singularity in Szekeres spacetimes can be 
"directionally naked" . 
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1 Introduction 



An important question regarding the final configurations resulting out of gravitational col- 
lapse of massive matter clouds is, in case a naked singularity develops rather than a black 
hole, whether it would be visible to a far away observer or not. This issue of local ver- 
sus global visibility of the singularity has important implications for the cosmic censorship 
principle. Whenever naked singularities formed, if they were only locally naked rather than 
being globally visible, that would mean the weak cosmic censorship would hold, even though 
the strong censorship may be violated [|I| . In its strong version the cosmic censorship states 
that singularities cannot be even locally visible, whereas the weak formulation allows for the 
singularities to be at least locally naked, but requires that they should not be globally naked. 

If the naked singularities formed but were always only locally naked and not visible to 
far away observers outside the event horizon, then though the strong form is violated, the 
validity of the weak version has the astrophysical implication that far away observers may 
never see the emissions from the singularity, or the very dense strong curvature regions 
near the same. This effectively seals the singularity inside the black hole. Thus, given a 
gravitational collapse scenario, it is of importance to determine if the singularity forming 
could be only locally naked, and if so under what conditions. On the other hand, if the 
singularity is globally visible, thus accessible to the outside observers, one would like to 
ascertain the conditions under which this would be possible, and to relate these to the 
nature of regular initial data from which the collapse develops. 

Our purpose here is to analyze a wide class of collapse models, namely the quasi-spherical 
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spacetimes given by Szekeres0, from such a perspective. Such an investigation is also of 
additional interest because these models are not spherically symmetric. The non-solvability 
of Einstein equations restrict most of the studies of formation and structure of singularities in 
gravitational collapse limited to spherically symmetric models, and a simple enough equation 
of state, namely p = 0. Though results on some specific solutions are available for more 
general equation of state p = kp^, still most of the work considering departures from 
spherical symmetry is numerical (see e.g. 0]). Analytic study of collapsing models with 
more general geometries is of considerable interest since strict spherical symmetry is a strong 
assumption. Thus we study here the formation of singularities in the Szekeres spacetimes, 
which do not admit any Killing vectors but allow an invariant family of two spheres, and 
has the Tolman-Bondi-Lemaitre (TBL) model as its special case. In a recent paper Joshi 
and Krolak |Q have shown existence of naked singularities in these models. Here we analyze 
the situation further to relate the formation of naked singularities or black holes with the 
initial data from which the collapse develops. The transition between the black hole and 
naked singularity phase is shown to be related to the nature of the initial data, and we also 
examine how the local versus the global visibility of the naked singularity changes with the 
change of essential parameters in the initial data. 

In section 2 we briefly review the Szekeres model. In sections 3 and 4, we derive the 
necessary and sufficient condition for the existence of outgoing radial null geodesies (RNGs) 
from the singularity, and we relate the formation of naked singularities with the regular data 
specified on the initial hypersurface. In section 5, we investigate numerically the conditions 
for the global nakedness of the singularity in the TBL and Szekeres spacetimes. The final 
section 6 outlines some conclusions. 



2 Quasi-spherical gravitational collapse 



The metric in a comoving coordinate system has the form, 

ds^ = df - M^dP - {dx^ + dy^), (1) 
The energy-momentum tensor is of the form of an irrotational dust, 

Qij = r^ij = g^i^i^ (2) 

where e{t, r) is the energy density of the cloud (here units are chosen such that c = SttG = 1). 
For convenience we adopt a pair of complex conjugate coordinates, ( = x + iy, ( = x — iy, 
in which the metric takes the form 

ds"^ = df - M^dr^ - N'^dCdC (3) 

Where M and N are general functions of t, r, ( and (. The integrated form of Einstein's 
equations is given by SzekeresP], 



Q(r,C,C) VT+7 ■ ' 



restricted to the case, iV' = ^ ^ 0, where the prime and dot denote partial derivatives 
with respect to r and t respectively. Here /(r) is an arbitrary function of r subject to the 
restriction / > — 1, and Q{r,(, C) is of the form 

Q = a(r)CC + B{r)C + B{r)C + c(r). (5) 



Here a and c are arbitrary real functions and S is a complex function of r, satisfying 



ac-BB^-, 5 = ±1, 0. (6) 



We will consider here the case S — +1. In this case it is possible to reduce the two-metric 
by a bilinear transformation of coordinates C) C to the form, 



R'{r,t)dCdC/\{CC + lf: (7) 



which, on subsequent introduction of polar coordinates ( — e^'^ cot ^9, ( — e~'^'^ cot ^9, 
reduces (7) to the normal form of metric on a two-sphere given by R^{r, t){d9'^ -\- sin^ 9dif^). 
Hence, the two-surfaces Sr,t defined hj t = const., r = const, are spheres of radius R{r,t). 
Here R{r, t) satisfies the "Priedmann equation" 



- fir) + (8) 



which is similar to the equation holding for the TBL models. Since we are studying collapse 
scenarios, we consider the collapsing branch of the solution, i.e. R < 0. One more free 
function, corresponding to the epoch of singularity formation, arises when one integrates the 
above equation and we get 

/ ^^^'^ + ^r'^'dR = t - toir) 

Functions F and / are arbitrary functions of r, and to{r) is another free function, which can 
be fixed using the scaling freedom. The density is given by. 



Therefore, the singularity curve is given by 



R{r,to)=0, {to>U), 



(10) 



Q{r, tsc, QR'ir, tsc) - R{r, tsc)Q'{r, tsc, C) = 0, (tsc > U)- 



(11) 



Here U, to and tsc give the times corresponding to that of the regular initial hypersurface, 
and the shell-focusing and shell- crossing singularities respectively. 

To analyze the shell-focusing singularity, it is essential to choose the initial data in such 
a way that ti < to < tsc- The range of coordinates is given by 



We can easily check that if the functions a(r), B{r), B{r), and c(r) are constants then 
the equation (9) reduces to the corresponding equation for the TBL metric. Amongst the 
functions a, c, B, B, F, f and to{r), equation (6) and a gauge freedom would leave total five 
free functions in general. In comparison to the TBL models, where we have only two free 
functions, namely F{r) (the mass function) and /(r) (the energy function), this model is 
clearly functionally more generic, allowing for a certain mode of non-sphericity to be included 
in the consideration. 



3 Gravitational collapse and singularity formation 



Ry^O, Q{r, t, C)R'{r, t) - i?(r, t)Q'{r, t, C) 7^ 



(12) 



Assuming the validity of Einstein's equations, the singularity theorems]^ ensure the for- 
mation of singularity in gravitational collapse if some reasonable conditions are satisfied by 



the initial data. For analyzing the detailed structure of the singularity curve, integrating 
equation (8) we get 



t - to{r) 



R^/^G{-fR/F) 

7f 



(13) 



where G{y) is a strictly convex, positive function having the range 1 >y > — oo and is given 

by 



f^-^, i>y>o. 



G{y) 



2 

3' 



y = o, 

^ ™VE^-^, 0>y>-oo. 



(14) 



Here to{r) is a constant of integration. Using the scahng freedom, we can choose 



i?(0,r)=r. 



(15) 



and this gives 



to{r) 



^^/''Gj-fr/F) 



(16) 



i.e. the singular epoch is uniquely specified in terms of the other free functions. The function 
io(r) gives the time at which the physical radius of the shell labeled by r = const, becomes 
zero, and hence the shell becomes singular. 



The quantity E! which is useful for the further analysis, can be expressed using equations 
(13-16) as 



R' = r"-i[(r? - I3)X + [e - (77 - ll3)X^'''G{-PX)] x [P + ^]V2] 
= r''-^H{X,r), 



(17) 



where we have used 



u = r", X = (R/u), r,(r) = r^, P(r) = r^, p(r) = r^, P(r) = pr'-\ 



j.a-1 py/2j.3{a-l)/2 ' ^3(a-l)/2 



;i8) 



The function (3{r) is defined to be zero when / = 0. AU the quantities defined here are 
such that while approaching the central singularity at r = they tend to finite limits, thus 
facilitating a clearer analysis of the singular point. The factor has been introduced for 
analyzing geodesies near the singularity. The constant a{> 1) is determined uniquely using 
the condition that Q{X, r) goes to a nonzero finite value in the limit r approaching along 
any X — constant direction. The angular dependence in the metric appears through ^, ^ and 
along radial null geodesies =const. 



Using the metric (3), the equation for the RNGs is, 



dt_ _ QR' - RQ' 



(19) 



Which in terms of R and u — r°' can be rewritten as 



dR 1 



du ar°' ^ 



. dr 



QR' \ oi ' 



U{X,u) + 



H{X,u) yJf + A/XRQ' 



a 



v^r+T QR'' 



(20) 



For spherically symmetric models we get Q' — 0, and the above equation reduces to the 



TBL case. The point m = 0, i? = is a singularity of the above first order differential 



equation. To study the characteristic curves of the above differential equation, we define 
X = R/r" = R/u. For the outgoing RNGs, dR/du > is a necessary condition. If future 
directed null geodesies meet the first singularity point in the past with a definite value of 
the tangent (X = Xq), then using equation (20) and I'Hospital's rule, we can write for that 
value Xq, 



Xn = lim — = lim — — = lim 



yr+T ^ QR" 



H{X,u) 



a 



(21) 



Since Q{0) ^ 0, Q'{0) = and near the central singularity R' oc R/r, if there is a definite 
tangent to the RNGs coming out, we can neglect [RQ') / {QR!) near the first singular point. 
Therefore the result is similar to the TBL modellHl, 



Xo= lim J7(X,r) = f/(Xo,0). 



(22) 



We have introduced here the notation that a subscript zero on any function of r denotes its 
value at r = 0. Defining 



V^(X) = f/(X,0) -X 



^ V^/o + Aq/X 



H{X,0) 



-X, 



a 



(23) 



the existence of a real positive value X = Xq for the above equation such that 



V{Xo) = 0, 



(24) 



gives a necessary condition for the singularity to be at least locally naked. 



In the neighborhood of the singularity we can write R = Xgr", where Xq is positive 
real root of above equation. To check whether the value Xq is actually realized along any 



outgoing singular geodesic, we consider the equation for RNGs in the form u — r°' — u{X). 
Prom equation (20) we have 



1 




1 ■ 




'—-X 




u 


. du 


u . 



U{X,u)-X + r\ J )^^^ 



(25) 



The third term on the right hand side goes to zero faster than the first near r — 0. The 
solution of the above equation gives the RNGs in the form u — u{X) . Now if the equation 
V{X) — has X — Xq as a, simple root, then we can write 



V{X) = {X-Xo){ho-l) + h{X), 



(26) 



where 



Aoi/o' 



XoNo 



is a constant and the function h{X) is chosen in such a way that 



(27) 



h{Xo) 



dh 



dX 



X=Xo 



0, 



(28) 



i.e. h{X) contains higher order terms in {X — Xq) and Hq — H{Xq, 0), Nq — N{Xo, 0) with 
N{X,r) = -R'r. 

As we are studying RNGs, Q and Q' are just the functions of r (i.e. u), therefore the 
equation for dX/du can be written as 



(Y Y\h^ ^ 
— (^A - Aoj — 

it II u u 



(29) 



in 



Here the function 



S = S{X, u) = U{X, u) - U{Xo) + h{X) + 



^f + A/X^XQ'u 



) 



(30) 



is such that S{Xq, 0) = 0. Note that unhke the TBL case, here S also has some functional 



those with fixed ( and (. Multiplying equation (29) by u and integrating we get 



where D is a, constant of integration that labels different geodesies. We see from this equation 
that the RNG given by D = always terminates at the singularity R — 0,u — 0, with 
X = Xq as the tangent. Also, for /iq > 1, a family of outgoing RNGs terminates at the 
singularity with X — Xq as the tangent. Therefore the existence of a real positive root to 
equation (24) is also a sufficient condition for the singularity to be at least locally naked. 

Near the singularity the geodesic equation can be written as 



Hence for Hq > 1 there is an infinite family of RNGs which terminate at the singularity in 
the past with Xq as a tangent. For ho < 1 there is only one RNG corresponding to D = 
which terminates at the singularity with Xq as a tangent. In {R,u) plane 



and for ho > 1 we have a family of infinitely many outgoing RNGs meeting the spacetime 
singularity in the past. For ho < 1 only D — geodesic terminates at the physical singularity 
(to(0), u — 0) in the past. 



dependence on ( and (. We can take care of this by considering radial geodesies, i.e. only 




(31) 



X-Xo = Du''°-\ 



R-XoU = Du^\ 



1 1 



A relevant question is, if only one RNG has the tangent X — Xq in the (u,R) plane 
then what is the behaviour of the RNGs which are inside this root and very close to the 
singularity? They cannot cross the root so they have to start from the singularity in the 
past. And if they are outside the apparent horizon, then certainly the only possible tangent 
they can have at the singularity is the same as the tangent made by the apparent horizon 
at the singularity. To check if our differential equation for RNGs can satisfy this near the 
singularity, we note that in equation (20) the term in the bracket goes to zero, while H blows 
up in this limit. In the limit to the singularity, we expect R/r^ to be Fq along the RNG. To 
the lowest order we assume that along the RNG 

and check for the consistency by determining Ms and 6 in the limit. For the f — case, 
using I'Hospital's rule we get. 



7-1 1- 

Fo = lim — = lim -—r- = hm 



F/R 



VX Fo 



(32) 



where X = R/r^. This leads to. 



5 — ?) — n and Ms — 



6^ 
nF„ 



That means that our assumption is consistent. One can also check that any other kind of 
behavior of the RNGs near the singularity is not possible. 
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4 End state of quasi-spherical dust collapse 



In this section, we analyze the end state of quasi-spherical dust collapse evolving from a 
given regular density and velocity profile. Our purpose here is to determine how a given 
initial profile influences and characterizes the outcome of the collapse. We will be mainly 
studying the marginally bound (/ = 0) case for convenience. The results can be generalized 
to non-marginally bound case using an analysis similar to that developed by Jhingan and 
Joship for the case of TBL models. 

The initial density profile for a fixed radial direction {( and ( are constant) is given by 
p = p{r), and if we are given the function Q then it is possible in principle to find out the 
function F{r) from the density equation (9). In order to give a physically clear meaning to 
our results, we assume all the functions (i.e. p, Q, and F ) to be expandable in r around 
the r = point on the initial hypersurface. This means that even though the functions are 
expandable, we do not necessarily require them to be C°° and smooth. Therefore we take 



oo 



P=Y1 Pnr 



(33) 



n=0 



oo 



Q = Y1 Q^^ 



(34) 



n=0 



where the regularity condition implies Qi = 0, and 



oo 



r 



,n+3 



(35) 



n 



n=0 



Therefore 



oo 



.n-1 



(36) 



n=2 



and, 



oo 



.n+2 



(37) 



n=0 



Substituting the values of these functions in the density equation and matching the coef- 
ficients of different powers of r on both the sides, we get the coefficients of function F{r) 
as. 



We note that Q has to obey the constraint given by equation (5) and (6) and F(r) does 
not have any angular dependence. This means that the angular dependence of different 
terms on right hand side of the equation (38) should exactly cancel each other. This leads 
to the conclusion that Pq, pi, p2 cannot have any angular dependence, and in general for 
n >S, pn to have some angular dependence at least p„_2 has to be non-zero. This confirms 
that local nakedness behaviour in Szekeres model is essentially similar to the TBL model. 

We now determine a and ©o which enter the equation for roots described in the previous 
section. Prom rj — rF'/ F we get rj as 




Pi 
4 ' 





(38) 



r]{r) ^3 + r]ir + r]2r^ + ri^r^ O(r^), 



(39) 
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where r^i, 7^2, ?73, ■■■■ are constants related to F{r). It turns out that we need the exphcit 
form of the terms up to order only. If all the derivatives p„ of the density vanish, for 
n < (q—l), and the q*'^ derivative is the first non- vanishing derivative, then T^, the g*'' term 
in the expansion for 77 is. 



" Fa 



(40) 



Here q takes the value 1, 2 or 3. In this case. 



(41) 



Prom equation (18), we get for marginally bound cloud {f — 0),l3 = p — P — and using rj 
from equation (41) we get. 



3 — 7) qFq 



3^3(a-l)/2 3^^^3(a-l)/2' 



(42) 



correct up to order q in r). The constant a is determined by the requirement that ©0 is finite 
and non-zero as r ^ 0, which gives. 



2q qFg 



(43) 



The root equation (23) with f — and rjo — 3, reduces to 



V{X) 



a 



1 - 




X + 



On 



-X, V{Xo)^0, 



(44) 



1R 



with 00 given by equation (43). The hmiting value of A = F/r" here is 



0, g < 3, 



Ao = < Fo, q = 3, 



(45) 



. oo, g > 3. 



As Ao takes different values for different choices of q, the nature of the roots depends upon 
the first nonvanishing derivative of the density at the center. 



This then leads to the conditions for local nakedness of the singularity in Szekeres models, 



similar to those obtained by Singh and Joshi 10 1 for the TBL models. The main similarity 



comes due to the Friedmann equation which is same in both the Szekeres and TBL cases. 
We consider the cases below: 



4.1 pi^O 



In this case, q=l, a = 5/3, and equation (44) gives. 



X, 



3/2 



Fi _3pi 
8po' 



2Fr 



(46) 



We assume the density to be decreasing outwards, pi < 0, hence Xq will be positive and 
singularity is at least locally naked. 



4.2 p,=0 , p^^O 



In this case, q = 2, a = 7/3, and equation (44) gives. 



X, 



3/2 _ F2 3^P2 

2Fo ~ lOpo' 



(47) 



Ifi 



Since density is decreasing outwards, i.e. p2 < 0, this implies Xq will be positive and again 
singularity is at least locally naked. 



4.3 pi = p2 = 0,p3 7^ 

Now we have g = 3, a = 3,Ao = Fq, Gq = —F^/Fq = — p3/2po, and pa < 0. As in this case 
Ao = -^0 7^ 0, the equation (44) for the roots becomes, 



^^1-M)(X- ) - X = 0. (48) 



Substituting X = Fqx^ and ^ = F^/Fq''^ = 3v^P3/2po^^ gives 

+ + - C = 0. (49) 

Now using the standard methods we can check for the possible existence of real positive 
roots for the above quartic. We see that for the values ^ < ^2 — —25.99038 real positive 
roots exist and we have naked singularity. In the case otherwise, we have collapse resulting 
into a black hole. 

4.4 = p2 = P3 = 

In this case of g > 4 we have a > 11/3, and Aq = oo. So positive values of Xq cannot 
satisfy the equation (44) for the roots, and the collapse always ends in a black hole (in 
the case of homogeneous collapse, q — oo and the collapse leads to a black hole, i.e. the 
Oppenheimer-Snyder result is a special case of our analysis). 
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4.5 a = 0, Ao 7^ 



In this case f = 0, F = rA, Aq 7^ 0, and the assumption in (34) is not vahd. Now from 
equations (17,18) we get, 



eo = | i/(x,o) = ^x + -^, (50) 



and equation (23) reduces to, 

x"^ + ^^qx^ ~x + ^0 = 0, (51) 
where X = x"^. 

We look for the real positive roots of this quartic in the same way as earlier. Substituting 
f3 = Ag^^/12, we get the condition for the real root, 

P> Pi = 17.3269 or p < p2 = 6.4126 x lO'l (52) 

But we can easily see that for the first case the above equation cannot be satisfied for any 
positive values of x. Therefore we get the condition for existence of naked singularity to be 
(3 < = 6.4126 X lO^"', which is the same as the condition obtained earlier by Joshi and 
Singh I 



IS 



5 Global visibility 



Amongst various versions of cosmic censorship available in the literature, the strong cosmic 
censorship hypothesis does not allow the singularities to be even locally naked. The spacetime 
must then be globally hyperbolic. Hence the case of existence of a real positive root to 
the equation (24) is a counter-example to strong cosmic censorship version. However, one 
may take the view that the singularities which are only locally naked may not be of much 
observational significance as they would not be visible to observers far away, and as such 
the spacetime outside the collapsing object may be asymptotically flat. On the other hand, 
globally naked singularities could have observational significance, as they are visible to an 
outside observer far away in the spacetime. Thus, one may formulate a weaker version of 
censorship, whereby one allows the singularities to be locally naked but rules out the global 
nakedness. In the present section, we now examine this issue of local versus the global 
visibility of the singularity for the aspherical models under considerations, and also for the 
TBL models, which provides some good insights into what is possible as regards the global 
visibility when we vary the initial data, as given by the initial density distributions of the 
cloud. 

In fact, all outgoing geodesies which reach r — Vc with R{t, Vc) > F{r^ will reach future 
null infinity. We assume /(r) and F{r) to be at least in the interval < r < Tc. Let us 
consider the situation when V{X) = admits a real simple root X — Xq, and a family of 
outgoing RNGs terminate in the past at the singularity with this value of tangent. Since the 
geodesies emerge from the singularity with tangent Xq, and the apparent horizon has the 
tangent Aq at the singularity, the local nakedness condition implies Xq > Aq. 

1Q 



Rewriting D in equation (31) in terms of Xc and Mc, where Xc and Uc are values of X 
and u on the boundary of the dust cloud, we have 



X — Xq — {X — Xc 



u 

Uc 



ho-1 



+ u''"-^ / Su-'"'du. (53) 



The event horizon is represented by the null geodesic which satisfies the condition Xc = A(rc). 
The quantity dR/du is positive for the outgoing geodesies from the singularity at r = 0. 
Therefore, all the null geodesies which reach the line r = Tc with Xc > A(rc), where the metric 
is matched with the exterior Schwarzschild metric, will escape to infinity and the others 
will become ingoing. The matching of the Szekeres models with the exterior Schwarzschild 



spacetime is given by Bonner [O . 



Our method here to investigate the global nakedness of the singularity will be to integrate 
the equation (20) for the RNGs numerically. For simplicity and clarity of presentation we will 
consider the case / = only. We expect that the overall behavior will be similar in / 7^ case 
also (for any given a), as evidenced by some test cases with / 7^ we examined; however, we 
will not go into those details here. The mass function F is taken to be expandable, and also 
we take the density to be decreasing outwards, as will be the case typically for any realistic 
star, which would have a maximum density at the center. There may be some possibility 
to study the behavior of RNGs near the singularity using analytic techniques, however, in 
a general case, we do not have analytical solutions for the concerned differential equation. 
Hence, we use numerical techniques to examine the global behavior of the RNGs, and to get 
a clearer picture, we keep only the first nonvanishing term after Fq in the expansion of F{r). 
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Thus, we have, 

F = For^ + F„r3+", 

i.e. for the TBL model, 

P = PO + Pn^-" 

where < Hereafter we fix Fq = 1 which is in some sense an overaU scahng. First 
we consider below the TBL model, and then some cases of Szekeres quasi-spherical collapse 
will be discussed to gain insight into the possible differences caused by the introduction of 
asphericity. 

5.1 TBL model 

Using the local nakedness conditions on the initial data as discussed in the earlier section, 
we analyze here the global behaviour of RNGs in the TBL model. 

5.1.1 Fi ^ 

In this case, the singularity is always locally naked, and we have a family of null geodesies 
coming out of the singularity with apparent horizon as the tangent at the singularity. We 
see in this case that the singularity can be locally or globally naked (Fig. 1(a), 1(b)). This 
basically depends upon the chosen value of pi, and where we place the boundary of the 
cloud. When the singularity is globally naked, all the RNGs which meet the boundary of 
the cloud with R < Rcru start from the singularity, and the other null geodesies start from 
the regular center r = before the singularity is formed. 

If we match the density smoothly to zero at the boundary, then Tc — — 3Fo/4Fi = — Po/pi, 

91 





Figure 1: The structure of singularity in the marginally bound (f=0) TBL model for pi 7^ 
case. Here po — 3, so the epoch of singularity formation is at t = to{0) — 2/3. (a) The 
event horizon meets the center r = below the singularity and no null trajectories escape 
to infinity. The singularity is only locally naked, (b) Singularity is globally naked. There is a 
family of null geodesies from the singularity which meet the boundary of cloud with > F^. 
(c) Graph showing transition from local nakedness to global nakedness, by varying pi. 
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and we see that for Fi < —1.12 i.e. pi < —4.48 the singularity is globally naked. In this 
case the event horizon starts at the singularity r = 0, t — 2/2) as shown in Fig. 1(c). It is 
seen that for the other case the event horizon formation starts before the epoch of formation 
of the singularity. It is thus seen that as we depart away from the homogeneous density 
distribution (the Oppenheimer- Snyder case), by introducing a small density perturbation 
by choosing a small nonzero value of pi, the singularity becomes naked. However, it will 
be only locally naked to begin with and no rays will go outside the event horizon, thus 
preserving the asymptotic predictability. It follows that even though the strong version of 
cosmic censorship is violated, the weak censorship is preserved. It is only when the density 
perturbation is large enough (Fig. 1(c)), beyond a certain critical value, then the singularity 
will become globally naked, thus becoming visible to the outside observers faraway in the 
spacetime. 



5.1.2 F2 7^ 

The overall behavior in this case (Fig. 2(a), 2(b)) is very similar to Fi 7^ case discussed 
above. The singularity is always locally naked whenever F2 is nonzero, and there is a family 
of RNGs coming out of singularity with apparent horizon as the tangent. If we match the 
density smoothly to zero at the boundary of the cloud, then Tc — y'— 3Fo/5F2 = \J —po/ P2- 
We see that for F2 < —3.1, i.e. p2 < —15.5, the singularity is always globally naked, while 
in the other cases it is not globally naked (fig. 2c). 

9"^ 





Figure 2: The structure of singularity in marginally bound (f=0) TBL model for F2 ^ 0, 

i.c P2 7^ case. We have po = 3, and = 2/3. (a) Singularity is only locally naked, 

(b) Singularity is globally naked. There is a family of geodesies from the singularity which 
meets the boundary of the cloud with R^. > F^. (c) The graph showing a transition from the 
locally naked to globally naked singularity by varying p2- 
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5.1.3 Fa 7^ 



This is the case when the first two derivatives of the density, pi and p2 vanish but pa ^ 0. 
In this case a = 3, and we have two real positive roots for the equation V{X) — if 
^ = F^/Fq^^"^ < —25.99038 = icrit- The singularity is then locally naked, and we see that 
it is also globally naked with a family of null geodesies meeting the singularity in the past 
with Xq as the tangent (Fig. 3(b)), where Xq is the smaller of the two roots. In the range 
— 1.5 > ^ > ^criti the apparent horizon forms at the time of singularity but the singularity 
is not even locally naked. Numerically we see that (Fig. 3(c)) in this case the event horizon 
formation starts before the singularity formation as expected. 

5.2 Szekeres models 

The spherically symmetric TBL models involve no angular dependence for the physical vari- 
ables such as densities etc. However, it may be possible to examine the effects of such a 
dependence by using the Szekeres models, which incorporate a certain mode of departure 
from the spherical symmetry. Thus, this may be a good test case to examine what is possible 
when non-sphericities are involved. With such a purpose in mind, we keep the form of the 
mass function F the same as that for the TBL models discussed above, and as a test case 
we choose the form of Q as given below, 

Q = cosh(ar)C( + j- — -. 

^ ^ 4cosh(ar) 

The constant a is chosen in such a way that all regularity conditions are satisfied. Again, we 
consider the subcases corresponding to different density distributions, as considered earlier. 
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Figure 3: Naked singularity in marginally bound (f=0) TBL model for F3 0, i.e. pa 7^ 
case, (a) Singularity is only locally naked, (b) Singularity is globally naked. There is a family 
of geodesies from the singularity which meets the boundary of the cloud with Rc > Fc. (c) 
The graph showing transition from globally naked singularity to a black hole by varying p^. 
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5.2.1 Fiy^O 



The singularity is always locally naked and there is a family of RNGs with R — F near 
the singularity. We choose a= 1.0. The general behavior of RNGs is similar to that of the 
similar case in TBL model. But there is a new feature emerging now, which is that the global 
visibility of singularity in some critical regions can depend upon the value of CC- That is, the 
singularity will be globally naked for a certain range of while for some other range of 
it would not be globally naked with the event horizon starting to form before the singularity 
for these values of CC (Fig- 4a). 

If we take the boundary of cloud at r = Tc = — 3Fo/4Fi then in the range — 1.06 < Fi < 
—1.01, the singularity is globally naked depending on the value of CC (Fig- 4a). Therefore 
the global visibility of singularity is directional. This feature distinguishes singularities in 
Szekeres spacetimes from TBL model in that the global visibility of the singularity will be 
related to the directional or angular dependence. For the values Fi < —1.06, the singularity 
is globally visible independent of the direction i.e. the values of C,C,, whereas for Fi > —1.01 
the singularity is not globally visible for any C,C, values. 



5.2.2 F2 7^ 

The singularity is always locally naked and there is a family of RNGs with tangent as the 
apparent horizon near the singularity. The general behavior of RNGs is similar to that 
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Figure 4: The graph showing the transition from locally naked to globally naked singularity 
in the marginally bound (f=0) Szekeres quasi-spherical models for different values of CC- 
Singularity is located at t = to(0) = 2/3. (a) Fi 7^ 0, a = 1.0, boundary of the cloud is at 
= -3Fo/4Fi. (b) F2 ^ 0, a = 2.0, boundary of the cloud is at Tc = (-3Fo/5F2)V2. (c) 
The globally naked singularity in the case F3 7^ 0, a = 3.9. 
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of Fi case above. If we take the boundary of the cloud at Tc = y-3Fo/5F2, then in the 
range —3.25 < F2 < —3.00 the singularity is globally naked depending on the value of 
(Fig. 4(b)). That is, there is a mode of directional global visibility for the singularity. For 
F2 < —3.25 the singularity is globally visible independent of the direction, and for F2 > —3.00 
the singularity is not globally visible for any values. 

5.2.3 F3 7^ 

In this case a — 3, and we have two real positive roots for the equation V{X) — if 
^ = Fs/Fq^^^ < —25.99038 = ^crit and the singularity is locally naked. We see in this 
case that if the singularity is locally naked then it is always globally naked. We can also 
see numerically that there is a family of geodesies meeting the singularity with Xq as the 
tangent (Fig. 4(c)). But still the equation of trajectories has some dependence on the angle. 

6 Conclusions 

We have examined here the global visibility and structure of the singularity forming in the 
gravitational collapse of a dust cloud, in the spherically symmetric TBL case, and also for 
the aspherical Szekeres models. An interesting feature that emerges, for the marginally 
bound case, is when we depart from the homogeneous black hole case by introducing a 
small amount of density perturbation which makes the density distribution inhomogeneous, 
even though the singularity becomes locally naked the global asymptotic predictability is 
still preserved. Such a singularity is not globally naked as our numerical results show, till 
the density perturbation increases and crosses a certain critical value, and thus for small 
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enough density perturbations to the original homogeneous density profile corresponding to 
the black hole case, the weak cosmic censorship is preserved. This, in a way establishes the 
stability of the Oppenheimer-Snyder black hole, when the mode of perturbation is that of 
the density fiuctuations in the initial data. It may be noted that this is a different sort of 
stability as opposed to the stability of the Schwarzschild black hole often referred to in the 
literature. What we have in mind here is the stability of the formation of black holes in 
gravitational collapse, involving the inhomogeneities of matter distribution, which remains 
invariant under small enough density perturbations in the initial density profile from which 
the collapse evolves. 

It remains an interesting question to be examined further as to whether this result will 
be still intact when we include the nonmarginally bound case, involving the nonzero values 
of the energy function /. 

Our analysis also confirms the earlier results [^], deducing in a more specific manner that 
under physically reasonable initial conditions naked singularities do develop in the Szekeres 
collapse space-times, which are not spherically symmetric, and admit no Killing vectors. This 
indicates the possibility that the earlier results for the final fate of spherical gravitational 
collapse might be valid when suitable generalizations to nonspherical spacetimes are made. 

As we have shown here, for the marginally bound TBL and Szekeres models, the global 
visibility of the singularity is related to the initial data. We classified the range of initial 
data for expandable mass functions F(r). We see that if the density decreases "fast enough", 
that is, when there is sufficient inhomogeneity, then the singularity is globally naked. In the 
Szekeres models, for some critical range of initial data we see that the singularity is globally 



naked only for a certain range of CC- Hence the singularity will be globally visible only in 
some directions. But there is no such directional dependence as far as the local visibility of 
the singularity is concerned. Thus our results point out and analyze the global nakedness of 
the singularity for the collapse of a dust cloud which is either spherical or quasi-spherical, 
with a regular initial data defined on a regular initial slice. 
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